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1 Introduction 

In this paper, we develop an approach to solve the statistical problem of an ideal gas in gen¬ 
eral 2-dimensional regions based on the conformal mapping method. Using the conformal 
mapping method and perturbation theory in Ref. [1], we calculate the spectrum of the ideal 
gas, and then we provide the partition function of the ideal gas on a general 2-dimensional 
region. 

It is difficult to solve the statistical problem within a general 2-dimensional region. 
The calculation of the partition function requires the spectrum of the ideal gas. In order to 
calculate the spectrum of the ideal gas in a general 2-dimensional region, we should consider 
a quantum billiard problem, which is defined as a region of space where a point particle 
is able to move freely, not being allowed to escape. We should calculate the eigenvalue 
of Schrodinger equation in this case. For a general 2-dimensional region, the calculation 
is very difficult. In fact, the exact solutions only exist for some special case, such as 
rectangular region, circular region [2], etc. We expect to develop an approach that can 
solve the statistical problem of the ideal gas with in a general 2-dimensional region based 
on some of these exact solutions. 
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There are some different approaches that have been developed to solve the Schrodinger 
equation within general 2-dimensional regions, you can get a review of this problem in the 
Ref. [2]. Among these approaches, the conformal mapping method is the most pronounced 
method, which is map the original region onto a given region by a conformal transformation 
[1], The boundary condition can be simplified by the conformal transformation, but the 
equation of motion become complicated. Further calculation of the eigenvalues requires the 
use of perturbation theory. So based on the approach elaborated in Ref. [1, 3, 4], we can 
calculate the eigenvalue on the general 2-dimensional region. Using the spectrum of the 
ideal gas, we can obtain the partition function of the ideal gas. 

The conformal mapping method is very important in mathematical physics. We know, 
that if the transformation is a holomorphic function defined on a region of the complex 
plane, then the transformation is a conformal map, which maps each point belonging to the 
region specihed onto another region. With the conformal mapping method, we can often 
simplify some physical problems. 

Now, we briefly review the conformal mapping method. Let us consider a Laplace 
equation over a 2-dimensional region D, 

(pxx T 4^yy — 0) (I'l) 

with Dirichlet boundary conditions. The points of the region D are denoted by z{x,y). 
According to the Riemann mapping theorem [5], we can map the region D onto a given 
region with simple boundary by a conformal transformation ui = f (z), i.e. 

u = u{x,y)+ iv{x,y) = f (z), (1.2) 


where z = x + iy. Under the conformal transformation, the equation transformed as 


(j^x T 4^uu T 2 (Ux'^x T ^y^y') 4^uv T {^x T '^y'} ^vv T i^xx T ^yy') 4^u T {Vxx T ^yy^ — 0. 

(1.3) 

Because the transformation to = f (z) is a, holomorphic function, u (x, y) and v (x, y) satisfy 
the Cauchy-Riemann condition, i.e. 


du _ ^ 

dx ’ 

du dll 

dy dx' 


(1.4) 


With Cauchy-Riemann condition, the equation on the new region can be simplihed. 


|/' (z)|^ {4>uu + (pvv) = 0. 


Obviously, this is still a Laplace equation, except for a conformal factor \f' {z)Y 


If we consider a Poisson equation 


(1.5) 


4^xx T (kyy — / {x, y) , (1-6) 

we can deal with this equation in the same manor as we did the Laplace equation. After a 
conformal transformation u = f (z), we obtain the equation on the new region, i.e. 

{u,v) ,y{u,v)]. 

1/ (^)l 


- 2 - 


(^UU T 4^vv 


(1.7) 



In Ref. [1-4], perturbation theory based on the conformal mapping method has been 
applied to solve the eigenvalue problem on general 2-dimensional regions. In this paper, 
we also calculate the partition function of ideal gas on general 2-dimensional regions by 
adopting the same approach. We will give numerical results of the partition function. 
Moreover the results of data fitting will be provided. 

The statistical thermodynamics of the ideal gas confined in a finite domain has been 
widely studied. In Ref. [6, 7], the authors calculate the size effect of ideal gases confined 
in a finite domain due to the wave character of atoms. For an ideal gas confined in a finite 
domain, there is a layer in which the density goes to zero which arises as a consequence arises 
as a consequence of the size effect. The boundary layer has been discussed in the Ref. [8, 9]. 
In Ref. [10, 11], the authors considered the effects of boundary and connectivity on ideal 
quantum gases in a confined space. In Ref. [12, 13], the authors discuss the properties of the 
hard-sphere gases in finite-size containers and convert the effect of interparticle hard-sphere 
interaction to a kind of boundary effect. In Ref. [14], the authors analyze the difference 
between boundary effects on Bosonic and Fermionic systems. Due to the size effect and 
boundary effect, we have to consider the non-extensive effect. The general discussion on 
the non-extensive effect can be found in [15]. In Ref. [16, 17], the authors discuss the 
mathematical basis of the boundary effects. Additionally, some related problems have also 
been studied. In Ref. [18], the authors discuss the thermal and potential conductivities 
by considering the size effects. In Ref. [19], the authors considered the Bose-Einstein 
condensation of the noninteracting bosons confined in a finite-size container. Based on the 
size effect, an ideal heat engine has been constructed in the Ref. [20-22]. 

In Sec. 2, based on the conformal mapping method and perturbation theory, we give 
the eigenvalue of a point particle which moving freely and not being allowed to escape from 
the boundary in general 2-dimensional regions. In Sec. 3, we provide the partition function 
of an ideal gas in a general 2-dimensional regions based on the eigenvalue spectrum we 
calculate in the Sec. 2. In Sec. 4, we shall give some examples, in which the numerical 
results and fitting results are given, and we compare the thermodynamic quantities of the 
different regions. The Conclusion is given in Sec. 5. 

2 Conformal mapping method and perturbation theory 

In this section, we briefly introduce how to calculate the eigenvalue of a point particle, that 
is moving freely and not being allowed to escape in general 2-dimensional regions. 

We consider a point particle freely moving in a general 2-dimensional region T> with 
complex boundary, namely a problem of 2-dimensional infinite potential well. But the 
region P is a small deformation of another 2-dimensional region D of which we have an 
exact solution. So we can deal with this problem by way of the conformal mapping method 
and perturbation theory in Ref. [2-4]. 

2.1 Conformal mapping 

We can describe a point particle moving in a region T> by using the Schrodinger equation 
with Dirichlet boundary condition. 
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( 2 . 1 ) 


2m 


■ 52 92 ■ 

du"^ dv"^ 


i>{u,v) 


Eil^ {u, v), 


For convenience, we set ^ = 1- We transform the region P onto region 11 by a holomorphic 
function, i.e. 


w = f{z), 


( 2 . 2 ) 


where w = u + iv, z = x + iy. According to the conformal mapping method we mentioned in 
the introduction, after the conformal mapping w = f (z), we obtain the equation of motion 
in the region H, i.e. 


- 


~ g2 Q2 - 

dx'^ 9y2 


'ip{x,y) = Ell; {x,y), 


(2.3) 


is the conformal factor, this equivalent local scaling factor after transfor- 


where S = 
mation. 

The boundary condition can be simplified under the transformation, but at the same 
time, the equation of motion is no longer represented by the Schrodinger equation. The 
Hamilton operator transforms as 


o = -s-i 


■ 52 52 ■ 

5x2 Qy2 




(2.4) 


where A = is Laplace operator. In order to continue our calculation of the 

spectrum, we require the use of perturbation theory. 


2.2 Perturbation theory 


In this section, based on the perturbation theory, we deal with the transformed equation 
(2.3). 

The operator O can be written as a symmetrized form [2], 


O 




(2.5) 


If the transformation we considered is a small deformation, then we can solve the problem 
by way of perturbation theory. 

For a small deformation, we rewrite the conformal factor as follows. 


S —)• = 1 + ya, (2-6) 

where a = T, — 1, y is a parameter denoting the order of the perturbation. Substituting Eq. 
(2.6) into Eq. (2.5), we expand the operator O in powers of y, i.e. 

O = Oq + yOi + y^02 + —, (2-7) 

where the coefficients of each power of y, 

Oo = -A, (2.8) 
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(2.9) 


Oi = [o-(-A) + (-A)o-], 

02 = 1 [2a (-A) a + Sa^ (- A) + 3 (-A) a^l , (2.10) 

O 

O 3 = [a^ (-A) a + a (-A) a^] - ^ [a^ (-A) + (-A) a^] , (2.11) 


Using perturbation theory, we may calculate the eigenvalue of operator O, the first two 
order read 

= Cn = ^nx,ny, (2.12) 

= - {n\ Oi\n) = -en^^ny {nx, %| a \nx, Uy) , (2.13) 

where \n) = \nx,ny) is the eigenfunction of zeroth order. We consider the correction to the 
eigenvalue up to first order, 

Enx,ny = ^nx,ny (1 “ {^x, ny \ a {Ux, ^j/)) • (2-14) 

Based on the conformal mapping method and perturbation theory in the Ref. [2], we 
obtain the perturbation expansion of the spectrum of a point particle moving in a general 
2-dimensional region. The following section we will deal with the statistical problem of an 
ideal gas in a general 2-dimensional region based on the result attained in this section. 


3 Statistical problem of ideal gas 


We already obtained the spectrum of a point particle moving in a general 2-dimensional 
regions in the last section by way of the conformal mapping method and perturbation 
theory. In this section, we consider the statistical problem of the ideal gas. 

Let us consider a system of ideal gas within a general 2-dimensional regions D. By 
definition we have the partition function [23], i.e. 

(3. 

Substituting Eq. (2.14) into Eq. (3.1), we obtain 


Z (fi) = = 




Z(/5) = EE exp [-I3en^^ny (1 - {nx,ny\a \nx, Uy))] . (3.2) 


Furthermore, we can calculate several thermodynamic quantities using the partition func¬ 
tion, such as internal energy U = — ^InZ (/3) [23]. 

Nevertheless, it is difficult to calculate the sum in the expression of partition function 
Eq. (3.2). We can adopt of numerical method to calculate the partition function in the 
specific case. The following, we will give some examples. 
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Figure 1. The square region 


4 Example 

In this section, we consider some specific examples. The examples considered in the fol¬ 
lowing are an exact solution— a solution with a square region, i.e. Fig. [1], The statistical 
problem of the 2-dimensional ideal gas in this region can be solved precisely. We have the 
spectrum on the square region, 


E-n — Cn — ^nx,ny — . {p^x T '^y) i 


(4.1) 


and the corresponding wave function is 


rirTT 


(x + 1) 


Sin 


(y + i) 


(4.2) 


We treat the result of the square region as the zeroth order when we deal with the following 
examples. We can obtain the partition function of ideal gas on this region, i.e. 




1 r 
4 


-!?3 0, e 4 




1 2 


- 1 


(4.3) 


where ds (0, z) is Jacobi theta function [24]. 

We can also calculate several thermodynamic quantities such as the internal energy 



















using the partition function, i.e.[23] 


C = -|lnZ(/,) 


^3 [0,e-^)-l 

The specific heat of the system is given by [23], i.e. 

d 


(4.4) 


Cv = ^U=^^U = 

dT 91 


- ) 9(5 


-U 


.ilA 




T2 -^3 


+ 


7?3 (^ 0 ,e 4 j - 1 

l?3 ("0,6-^^ -1 


(0,1) / „ -2-^ 


0, e 


+ e 4 


i 93 0,e-— -1 


1?. 


( 0 . 1 ) 


.IZA 


I 0 , e 4 


0 , e 4 

(4.5) 


In the sections that follow we shall consider various cases. 

4.1 Transformation of w = f (z) = z + 

Let us consider the 2-dimensional regions Di with complex boundary. Fig. [2] where we set 



A = 0.1. After conformal transformation, we map the regions Pi onto regions P Fig. [1]. 
4.1.1 Spectrum of eigenvalue 

First, we calculate the conformal factor of the transformation oi w = f (z) = z + Az^, i.e. 

2 


S = 


dz 


= 1-1- 4xA -|- 4A^ (x^ -|- y^) , 


(4.6) 
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then we have 


(T = S — 1 

= 4Ax + 4A^ {x^ + y^) 


(4.7) 


Substituting cr, zero order eigenvalue Eq. (4.1) and zero order eigenfunction Eq. (4.2) into 
the expression of first order correction to eigenvalue Eq. (2.13), we obtain 
.2 /-l /•! 


^ = -^{nl + nl) 

= + (l- 


' dx J dy sin^ 
3 


najTT 


TT^n^ 


L 2 

3 

vr^riy 


(x + 1) 


sm 


L 2 


{y + 1) [4Ax + 4A^ (x^ + y^)] 


(4.8) 


4.1.2 Partition function 


Now, we calculate the partition function of 2-dimensional ideal gas on the region Pi. Sub¬ 
stituting the eigenvalue into Eq. (3.2), we have 


€ 


1 - -A^ 1 - 


TT^n^ Tr'^ri^yj 


(4.9) 


Obviously, it is hard to give an analytical result of the sum, so we can calculate the numerical 
result of the sum for different (3 (namely different temperature). We calculate the partition 
function by choosing remove different values of /3. We list numerical results in the Tab. [1], 


Table 1. Numerical results 


/3 

0.01 

0.02 

0.03 

0.04 

0.05 

0.06 

0.07 

zm 

27.158 

12.5049 

7.80701 

5.5305 

4.20164 

3.33767 

2.73481 


0.08 

0.09 

0.10 

0.20 

0.30 

0.40 

0.50 

Z{/3) 

2.29255 

1.95574 

1.69171 

0.594704 

0.288229 

0.15799 

0.0917356 

/3 

0.60 

0.70 

0.80 

0.90 

1.00 



Z{/3) 

0.0547901 

0.0331811 

0.0202314 

0.0123765 

0.00758345 




In principle, we can give the numerical result for different temperature. But we still 
give an expression of the partition function based on the table of numerical results. The 
perturbation calculation is based on the square region P, so we can set the partition function 
on the region Pi as same as the form of the partition function on the square region H, 

1 2 


^(/3) = (0,6-?-^/^)-!' 

Taking the values from the Tab. [1], we give the htting value of (m,n), as 


(4.10) 


m —^ 1.0317, n ^ 1.004 


- 8 - 





























so we obtain the partition function of 2-dimensional ideal gas on the region Pi 


zm = 


1.0317 


. / 1.004 _2 

0, e 4 


7r2/3 ) _ ^ 


1 2 


(4.11) 


The internal energy of this system is given by [23], i.e. 


d 


C = -^lnZ(« 

O / 1.004 _2 

l.OOdTT^e ’ (3,0,6 ^ 


2(^93 (0,e-^-^^ 

We also have the specific heat of the system [23], i.e. 

d 


- 1 


Cv = = —^U = 


r) 1.004 ^2 Q A _o 

1.004^6 2 ^ f^TT^/3^ 


u 




(o,e 


1 - 1 


(4.12) 




+6 4 


i7r2/3 


193 ( 0,6 4 


-1 




+ 


)}■ 


^ / 1.004 _2o\ 

193 (0,6 4 1 


( 0 , 2 ) 


i9. 


(4.13) 


0,6 4 


The effect of the thermodynamic quantities after the transformation are shown in the 

Fig- [3], 



(a) The internal energy (b) The specific heat 

Figure 3. The solid line represent the thermodynamic quantities of the square region, and the 
dotted line represent the thermodynamic quantities of the Pi. 


4.2 Transformation w = f {z) = z + \e^ 

Let us consider the 2-dimensional regions T >2 with complex boundary Fig. [4], where we set 
A = 0.1. After conformal transformation, we map the regions P 2 onto regions 11 Fig. [1]. 
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4.2.1 Spectrum of eigenvalue 

First, we calculate the conformal factor of transformation oi w = f {z) = z + Xe^, i.e. 

2 

= 1 + 2e^Acosy + (4-14) 

then we have 



fj = S — 1 

= 2e*Acosy + A^e^*. (4-15) 


Substituting a, zero order eigenvalue Eq. (4.1) and zero order eigenfunction Eq. (4.2) into 
the expression of first order correction to eigenvalue Eq. (2.13), we obtain 

=-^ [nl + Uy) j dx J dysin^ (x + 1) sin^ (y + 1) (2Ae'^ cos y + 


= -^{nl + nl) 


(e — e vr^n^n?, sin 1 „ 

^ ^ a —TV + A^ 


L 2 

TT^n^ sinh 2 


(vr^n^ + 1) (vr^nl — l) 2 (vr^n^ + 4) 

4.2.2 Partition function 


(4.16) 


Now, we calculate the partition function of 2-dimensional ideal gas on the region ^ 2 - Sub¬ 
stituting the eigenvalue into Eq. (3.2), we have 


Z(/3) = ^e-^^ 

e 





/ _ (e-e TT^n^ng sin 1 
V (7r2n2 -h 1) (7r2n2 - l) 


, 9 Tr'^n'i sinh 2 \ 

(4.17) 
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Sililar to case 1, we calculate the numerical result of the sum for different /3 and give the 
numerical results in the Tab. [2], 


Table 2. Numerical results 


/9 

0.01 

0.02 

0.03 

0.04 

0.05 

0.06 

0.07 

zm 

34.458 

16.0312 

10.0942 

7.20541 

5.51287 

4.40864 

3.63556 


0.08 

0.09 

0.10 

0.20 

0.30 

0.40 

0.50 

Z{P) 

3.06661 

2.63197 

2.29021 

0.852253 

0.437876 

0.255918 

0.15969 


0.60 

0.70 

0.80 

0.90 

1.00 



Zif3) 

0.103308 

0.0682117 

0.0455631 

0.0306348 

0.020674 




Taking the values from the Tab. [2], we obtain the partition function of 2-dimensional ideal 
gas on the region 7)2, 


z (13) = 15^ [tfs (o, - l] t 


The internal energy of this system is given by [23], i.e. 


d 




. / 0.788871 _2/3\ 

0,e--1 


We also have the specific heat of the system [23], i.e. 

d 


Cv = ^U = ^^U = 

0.788871 2 


<9/3 


-U 


O U.y888yi_2/3 A 

_ 0.788871^6-f (0,1) 


8 

c_ 

H-e 4 


I^s (^0, e 4 

0.788871 ^2 ^ 


0.788871 


- 1 


(4.18) 


(4.19) 


_ 0.788871 ..,2 

193 0, e 4 


^(0.1) _ jJ ,,(0,2) 

)}. (4.20) 


0.788871 2, 


^ / 31 _ 1 


( 0 , 1 ) 02788 ^,^ 2/3 > 


(0,e-^ 


The effect of the thermodynamic quantities after the transformation are shown in the Fig. [5], 


4.3 Transformation o^ w = f {z) = z + \ sin z 

Finally, we consider the 2-dimensional regions T>^ with complex boundary Fig. [6] , where 
we set A = 0.1. After conformal transformation, we map the regions T>^ onto regions fl Fig. 
[ 11 - 
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(a) The internal energy (b) The specific heat 

Figure 5. The solid line represent the thermodynamic quantities of the square region, and the 
dotted line represent the thermodynamic quantities of the 2 ? 2 - 


1.5 



l.O 

e.5 



1.5 

-r 

o 

1 

<> 

0.5 I.C 

l.i 


-0.5 




-l.O 


. 

• 1.5 




Figure 6 . A = 0.1 


4.3.1 Spectrum of eigenvalue 


First, we calculate the conformal factor of the transformation of w = f (z) 

i.e. 


S = 



1 + 2A cos X cosh y + (cos 2x + cosh 2y ), 


z + Asinz, 


(4.21) 


then we have 


fj = S — 1 


= 2A cos X cosh y + - A^ (cos 2x + cosh 2y). 


(4.22) 
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Substituting a, zero order eigenvalue Eq. (4.1) and zero order eigenfunction Eq. (4.2) into 
the expression of first order correction to eigenvalue Eq. (2.13), we obtain 
^2 


^ + nl) j ^dx j ^ 


dy sin^ 


UrTT 


{x + 1 ) 


sin 


(y + i) 


2A cos X cosh y + -}? (cos 2x + cosh 2y) 




(e-e vr^n^n^ sin 1 2 '^^’^xSin 2 „'K‘^n^s\ii\i2 

(vr^n^ — 1 ) (vr^n^ + l) dvr^n^ — 16 dvr^n^ + 16 


(4.23) 


4.3.2 Partition function 


Now, we calculate the partition function of 2-dimensional ideal gas on the region P 3 . Sub¬ 
stituting the eigenvalue into Eq. (3.2), we have 


Z(/l) = ^e-^^ 

(e - e“^) TT^n^n^sinl 2 ^^^sSin 2 ^'n‘^Ti?ySmh.2 
(vr^n^ — 1 ) (TT^nl -|- l) dvr^n^ — 16 dvr^n^ -|- 16 

(4.24) 

Same as case 1 , we calculate the numerical result of the sum for different (5 and give the 
numerical results in the Tab. [3], 



Table 3. Numerical results 


/3 

0.01 

0.02 

0.03 

0.04 

0.05 

0.06 

0.07 

zifd) 

34.1644 

15.8955 

10.0092 

7.14507 

5.46694 

4.37208 

3.60555 

13 

0.08 

0.09 

0.10 

0.20 

0.30 

0.40 

0.50 

Z{f3) 

3.0414 

2.61042 

2.27152 

0.845474 

0.434401 

0.25384 

0.158328 

P 

0.60 

0.70 

0.80 

0.90 

1.00 



Z{(3) 

0.102364 

0.0675327 

0.0450658 

0.0302677 

0.0204026 




Taking the values from the Tab. [3], we obtain the partition function of 2-dimensional ideal 
gas on the region P 3 , 


Z(/3) = 


0.996072 


■da 0 , e 4 


0.787783 ^2^ 


- 1 


The internal energy of this system is given by [23], i.e. 

C = -|lnZ(« 

0.787783vr2e-^-^/^#^) ^3,0, 


-!?3 0 , e 4 


0.787783 


- 1 


(4.25) 


(4.26) 


13 - 











































We also have the specific heat of the system [23], i.e. 


d 


d 


Cv = ^U = -^U = 

dT Q 


d 


d/3 


-U 


0.787783% 


n 0.787783 _2o ..n 

2p-5- tt 


_ 0.787783 2. 

0,e 4 ^ I- 


1 2 


- 1 


-7? 


( 0 , 1 ) 


0 ,e 


0.787783 2/ 


+ 


0.787783 2 


H-e 4 


77^/3 


. / 0.787783 _2o\ 

7?3 (0,e- n - 1 


& 


( 0 , 1 ) 


0 , e 


0.787783 2 


7r'^/3 


i?3 (^0,e 4 

)}■ 


0.787783 2 


- 1 


7? 


( 0 , 2 ) 


(4.27) 


The effect of the thermodynamic quantities after the transformation are shown in the 
Fig. [7], 



(a) The internal energy (b) The specific heat 

Figure 7. The solid line represent the thermodynamic quantities of the square region, and the 
dotted line represent the thermodynamic quantities of the 'D^. 


5 Conclusion 

In this paper, based on the conformal mapping method and perturbation theory, we present 
an approach to solve the statistical problem of ideal gas within general 2-dimensional re¬ 
gions. We considered some example and give numerical results. For more general case, we 
can state the transformation as follow: 

f{z) = z + '^Xnz'^ -h -h ^/3„sinnz. (5.1) 

Then we obtain the transformation by calculating the undetermined parameters A^, an and 
/3n- Using the approach given in this paper, we are able to solve the statistical problem of 
ideal gas in general 2-dimensional regions. 
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